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Математичне моделювання перистальтичних процесв 
на клиткових нейронних мережах 


В статье рассмотрены методы и средства моделирования перистальтических процессов в биологических 
объектах. Постановка задачи представлена системой, состоящей из уравнения Навье-Стокса и уравнения 
непрерывности потока. Для решения задачи предложено использовать локально-асинхронный метод, 
ориентированный на решение уравнений математической физики в клеточных нейронных сетях. 
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Введение 


Перистальтические процессы, или процессы продвижения жидкости в трубчатых 
поверхностях за счет колебания стенок этих поверхностей, распространены в биологи- 
ческих объектах. В организме человека, например, они играют важную роль в сердечно- 
сосудистой и пищеварительной системах [1]. Создание математических моделей пери- 
стальтических процессов дает возможность изучить физическую природу этого явления 
с целью компенсации и нивелирования негативного влияния на живой организм, 
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вызванного отклонениями параметров упомянутых процессов. Математическая модель 
представлена системой уравнений, состоящей из уравнения Навье-Стокса и уравнения 
непрерывности потока [2]. Решение данной системы в аналитическом виде не пред- 
ставляется возможным, а использование стандартных численных методов весьма затруд- 
нено при рассмотрении соответствующей краевой задачи в нестационарном варианте 
в трехмерной области с конечными геометрическими размерами и переменными 
формами. В то же время математическая модель в такой постановке вызывает значи- 
тельный практический интерес, что делает актуальными исследования по созданию 
новых методов и средств, позволяющих получать результаты исследования модели в 
реальном масштабе времени. 

Целью исследований является повышение эффективности моделирования пе- 
ристальтических процессов путем развития параллельных асинхронных методов, ори- 
ентированных на реализацию в параллельной вычислительной среде со свойствами 
клеточных нейронных сетей. 


Постановка задачи 


В рамках математической модели перистальтического процесса поток жидкости 
в трубчатом объекте будем рассматривать как ламинарный поток вязкой несжи- 


маемой жидкости, ограниченный во времени отрезком Ёе 0,7 ] и в пространстве 
областью Ос 83. Основные характеристики данного потока представлены полем 


скоростей и:Ох 0,7 |-> В? и давлением р:©@х 0, Г ] —В, для которых взаимная 
зависимость и зависимость от граничных условий определяется системой уравнений, 
состоящей из уравнения Навье-Стокса и уравнения непрерывности потока. Рассмотрим 


систему декартовых координат (х, 9,2). Тогда й=м+ №+ми Е=Ш х+ЛВу АЕ. , 


где (1, У.) — единичные векторы, параллельные соответствующим осям координат. 


Используя эти обозначения, представим систему уравнений в координатной форме [3]: 


ди 1 бр _п 92и д?и , д’ „9 ди ди 
д рах р| а? ду? 2022 0х д 62 


д» Тр п| 02» 02» 0] № 65 @ 
Е = + + и У и ЧА 
и рб Р| 0х? ду? 02} хх д 4 а 
ди» Тор п 02» 02» 02» би 0 0 
= яр и У ЕР 
и р р| д ду? а=2 0х ду 02 
0% _| 
х д & `` 
р- плотность, ий — динамическая вязкость, Ё — поле объемных сил. 
Определим начальные значения скоростей при #=0 
и=ио(х, у,2)= 0; у= 0 (х, у,2)= 0;ж= мо(х, у,2)= 0. (2) 
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Учитывая особенности перистальтических процессов в биологических объектах, 
при исследовании этих процессов будем использовать граничные условия, опреде- 
ляемые скоростью движения ограничительной поверхности. Они задают значения 
искомой переменной в непосредственной близости к поверхности при отсутствии 
скольжения вдоль поверхности. Стенки поверхности формы движутся со скоростью, 
направление вектора которой совпадает с нормалью и к оси 2, а величина модуля 
О и определяется функцией: 


ОИ (ь2,6,ф,4,7)= О тах О л)-И 2(2,6,6,7), (3) 


где Итах - величина максимально допустимого значения модуля скорости, 
Я — величина временного периода изменения скорости деформации, д — порядковый 
номер пространственного шага, у -— величина шага смещения деформации вдоль оси, 
ф- период деформации вдоль оси. 

Временную составляющую функции изменения модуля скорости будем опре- 
делять из выражения: 


были? и! 4) 


Два периода функции И’; задают полный цикл колебания поверхности. Первый 


период соответствует росту деформации ограничительной поверхности трубчатого 
объекта, а второй период описывает процесс возвращения деформированной поверх- 
ности до первоначального состояния. 


График временной составляющей функции изменения модуля скорости И’, по- 


казан на рис. 1. 


и») 


Рисунок 1 — Временная составляющая изменения скорости деформации 
поверхности в перистальтическом процессе 


Пространственный шаг изменяется во времени в соответствии с выражением 


Выражение для пространственной составляющей И 2(2, ‚ф,у) задает форму пе- 


ристальтических колебаний. При перемещении жидкости за счет продвижения вдоль 
оси 2 зон сжатия и расширения, пространственная составляющая И Е: ‚ф,у) при- 


обретает вид: 
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< 9=0, 


— 1” вии? [1-6] прио<2<8 
д 


О. (2,6,$,у)=4-зшт [Е (упр 88 О<б<д а, (5) 


| | ЗФ 

7 [АЕ-ви вой [==- 6-1 р —<2<29 

— 517 [Ев И. 
ф 2 


—5т РЕВ р а = 
ф 2 7, 


5? [ев при * < 2<29 


График изменений пространственной составляющей И „ ‚ представленной выра- 
жением (5), показан на рис. 2. 


02 (2, 5, Фф, у) 


Рисунок 2 — Пространственная составляющая изменения скорости деформации 
поверхности в перистальтическом процессе типа «сжатия и расширения» 
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Дискретизация краевой задачи 
Разностная схема для и(т + 1) : 


2 2 2 
и(т+1=и()+& ие ие и т та и Ее те : (6) 
Р\ 0х 0’ 02 


Введем замену: 


т 0?и ди 02и „ди _ ви _ „9 
0х ду 02 


у и “в (7) 


Учитывая, что значение А(т) связано с временным уровнем т, из (6) получим 
соотношение: 


и(г +1)= 4-Я 2+1) (8) 
р’. . © 


Выполним преобразования (6) и (7) для остальных координатных уравнений и 
получим выражения для у(т+1) и и Не 1): 


№+-в6)-9. 2+9, Дт +0- сб) Я ЕН. (9) 


Подставив выражения (8) и (9) в уравнения непрерывности потока, получим 
уравнение Пуассона для давления: 


97 р( +1) р 97 р( +1) у 92 ре +1) _р дА(г) й ОВ(т) : ОС(т) 
ах? ду? аг2 &\ дх ду 02 }` 


Дискретизируем левую часть уравнения Пуассона (10) и подставим дискретные 


(10) 


значения 4; ‚к, В; к›Ср;к в его правую часть. Как результат получим систему урав- 


нений: 
рык @+-2р; ке +1)+ рак +1) Е 


(5)? 
Рак @+-2р:(е+1)+ р; 1+1) ь 
(5 
Рид @+Ю-2р; ке +1)+ рр +1) й 
(=) . 
= 2 (т) - А, (т-1) и Вых (т)- В (7—1) к Си (т)- С, (т я ’ 0 
01 ох ду 02 


где 1=1.--*,Йпах» /=Ъ-"*, Лпах › А =Ъ**",Ютах . 


Дискретные значения граничных условий для р, А, ВТС: 


Ро, 7.К = Ру уК› Р: + 7,.К = Р: „1,К? 
| Л Л тах] тпах >. ее р (12) 


Ао, = и, „К 2 Чи = ЕЙ 2 
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Рок =РАК› Рук = Ррриах,К> 
Лпах Лтах р ах ее (13) 
В;,0,к > и, 0,К 2 ВНЕ =: ах „К 7 


Р+,},0 | РЫ Ру, ушах Е. Ру, Ктах й В - : с 
(@ = С = =1,*‘,Йпах, Л=Ъ""‘,Лтах. (14) 
1.130 — 91, 1.05 ах. Ча? 


Граничные условия влияют на вид системы уравнений (11) при использовании 
ее для определения давления в центрах элементарных кубов @ трехмерной сетки, 
расположенных на границе исследуемой области. Система уравнений (11), учитываю- 
щая упомянутые модификации, имеет вид: 


И (р, ИР Е +1)+ т (1 Ее - рык (Е +1). 


(55)? 
9” (>. ны +1- рык +1) ати (>. ры (+1- рр +1). 
(5 
в“ (>, Це +- рык +1)+ о ы (>, е+- рык +1)) 9 
= = (5) 
_р аа -9. ВоВ -9. Се — 
р бх бу 5 


где #=1,--*„Йнах › Л=\-**, Ллпах › К=Ъ--*,Ющах . 


шт т Г=Ь пах -| Г < Плах› 


[72% а ре г={1, К}. 


ОЕ 


Система линейных алгебраических уравнений (15) состоит из йлах ЛлахАтах 


уравнений с неизвестными, составляющими множество с мощностью йдах/тахАтах : 


р 


ыы] а ыы 2 


Следовательно, такая система может быть решена с использованием одного из 
итерационных методов. Выбор такого метода является нетривиальной задачей, учиты- 
вающей также и особенности средств, планируемых к применению для ее решения. 


Локально-асинхронный метод 


Представим систему уравнений (15) в виде: 
Гу(х}=Е(х), (16) 


ее 


Будем решать краевую задачу методом, основанным на принципах асинхрон- 
ности и локальности [4]. Главная особенность данного метода заключается в том, что 
он предполагает реализацию вычислений в параллельных структурах, обладающих 
свойствами клеточных нейронных сетей [5]. 
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Клеточную нейронную сеть, ориентированную на данный метод [6], можно пред- 
ставить как конечную совокупность клеточных нейронов, реализующих глобальный 


[©.6) 
оператор Г. Асинхронная итерационная последовательность [7] {Р(т)} ‚1 Векторов 


р(т) = [р ТИ (тах оах мах Е В формируется в соответствии с выражением: 


рр (т), (5 1,Ю Ум, 


1 2 
РК | К (")} Ру, к [2 1) 


Рьлк(т+1)= 3 4 3 лы о 
м Рори лет рак ЕЕ (т) |, ЛЮ, 
5 6 3 
где $= | р (т) 1=1,...,6; т= 2... — множество задержек передачи данных 


в нейрон с координатами (г, /,К) от окружающих его нейронов на итерационном 
шаге т, р(0)= [р т (ооо для Атях — вектор, определяющий исходное состояние 
Эа — 


распределения давления. 
Дополнительные условия, задающие свойства вычислительного процесса [8]: 


5= шах (м-5, „)) И О Ы 
Бу, (Ь1,К) а 
5ЕМ№д; 


а БО о оных 198 МОУ ЕМ, 
Разностная схема для локально-асинхронного метода (17) имеет вид: 


И [ок | 


о Пр, вЫ ик) др кб? к) ат р, ры ик) мер а к) 


(5 (5, 


тт 


С, РЕ в (т)) ча, Вены в (т)) 


го —_ 4.0 1 2 3 4 > 6 — 
Где 5; (+) = о. А ЗА множество задержек 


передачи данных от элементарных кубов: 
|, >, >, К >91, 1-Е, + › 9, к-р, „К } ; 
Лк - функция правой части; 


у — релаксационный параметр, находящийся в интервале 0,2] . Для данной крае- 
вой задачи принято значение у =1,7. 


т =1,...-— текущий номер итерации на временном шаге т. 
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Для решения рассмотренной краевой задачи будем использовать трехмерные кле- 
точные сети со структурой послойного циркулянтного графа. Простейшая из таких 


структур содержит 64 нейрона и представлена графом С(64;1,4,16) на рис. 3. Такая 


структура недостаточна для практической реализации упомянутой математической 
модели, однако дает представление о характере межнейронных связей и взаимном 


пространственном расположении нейронов. 


ПИ 


о 
а с ЧА т 
КМК 
к Е В 
Чо № 
Е 
тео 
НЕМ АНИ Ь 
тыщ 
к. а ПУ ТЫ 
АИ еМ 
КРИ РЕ р 
В ЕАК 
М Е ЕЕК А 
к Ц 
Ау НЫ АВК 
Е 
АЕ 
К 1 т МАЕ ВА 1 
Е 
ИХ 
вом а ЕВ 
Е АКА 
Ш КЕМ 
| А 
| © Е а 
^^ Е 
х ОУ 


Рисунок 3 — Граф С(64;1,4,16) трехмерной клеточной нейронной сети 


Связи между нейронами клеточной нейронной сети вдоль осей х, у и 2 обра- 


зуют совокупности вершин, которые могуг быть достигнуты путем выполнения шагов: 


(х=1тоа1б, (у=4)тод16 и (2+16)тоа 64. 


Существует строгое доказательство возможности существования подобных струк- 
тур произвольного размера при условии, что количество нейронов, расположенных 


вдоль одного ребра, равно 2" ‚где и=2,3,4,---. 


Клеточный нейрон представляет собой специализированное вычислительное уст- 
ройство, обеспечивающее выполнение разностной схемы (18) и обмен информации с 
соседними нейронами в соответствии с локально-асинхронным методом (17). Прин- 
ципиальным отличием клеточного нейрона от специализированного вычислительного 
устройства является присущее ему свойство обучения, обеспечивающее увеличение 
скорости сходимости итерационной последовательности при решении краевой задачи. 
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Формирование итерационной последовательности в соответствии с локально- 
асинхронным методом (17) допускает произвольное сочетание данных, полученных от 
соседних нейронов. Такой подход фактически эквивалентен использованию на каждом 
итерационном шаге нового численного метода из числа известных синхронных числен- 
ных методов. Очевидно, что скорость сходимости в данной вычислительной ситуации 
для каждого такого метода будет различной. Таким образом, суть обучения клеточного 
нейрона состоит в запоминании для каждого конкретного решения краевой задачи 
такой итерационной последовательности, которая обеспечивает максимальную ско- 
рость сходимости локально-асинхронного метода в целом. 


Выводы 


В данной работе рассмотрен подход к моделированию перистальтических про- 
цессов в трехмерной постановке на основе решения системы уравнений, состоящей 
из уравнения Навье-Стокса и уравнения непрерывности потока. Поскольку решение 
такой системы при использовании сложных движущихся границ весьма затруднено, 
предложен комплекс методов и средств, позволяющих получать результаты моделиро- 
вания в реальном масштабе времени. В качестве средства реализации математической 
модели использована клеточная нейронная сеть со структурой послойного циркулянт- 
ного графа. Свойство обучения клеточного нейрона позволило обеспечить повышение 
скорости сходимости локально-асинхронного численного метода. 
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Математическое моделирование перистальтических процессов... 


ТН 


КЕЗОМЕ 


М.А. МоуоагуМу, В.В. МеяегепКо 
эпишайоп о{ Релящис Ргосеузех Иятя 
Се|шаг Меига Ме лотку 


ТЬ15 рарег 15 деуоде4 №0 деуе]ортепЕ оЁ Фе тео4$ ап4 10015 Юг зпишайоп рег1з- 
{ас ргосеззез ш БМоог1са| обес. Зиич]айоп 15 регюгте4 Бу оМашше а зедиепсе 
ОЕ зоНоп$. Еогианоп оЁ Ше ргоет ргезещеа Ъу аезужет оЁ Мауег-ЗюКез ап4 
сопипийу едааноп5. Оие ю Фе Ёс! {а Фе зпашайоп 5е5 а зорзисаеа Югт оГа гег1оп 
У тоуше Боипдатез, й 15 уегу АЁЯсий ю изе га опа| тео ап4 {0015 {0 сеё е 
гези $ оРтезеагсН т геа]-@ите. 

п огдег о зрее4 пр гезаЁ оепеганоп оЁ зпиайопз, се|шаг пепга! пебмок ууй @гее- 
Аитепз1опа| 1ауеге4 зёгасваге оЁ сисшапЕ огарН 1$ зиооезе4. СеПааг пеигоп$ аге юсаёе4 
а Бе 5Иез оЁ 1$ сеШ\аг пейга] пебмогК. ТВезе аге зреслаПте4 сотрийпе 4еу1сез УВ 
фе ргорему оЁ]еагише. ТБе Баз1с сотршайопа| тео4 15 е |осаПу азупсгопои$ тео4. 
ТЬ$ те@о4 1$ изе4 10 рго\4е Ше азупс№гопои$ пааге оР соттишсаНоп ргорег © фе 
се|аг пепга! пебмотКз. Ргорегиез оЁ{1е |осаПу азупсгопои$ тео4 а150 ипаегИе се|аг 
пеигоп$ наше. ТВе аиогз Вауе изе4 {1е ргорегие$ оЁ {11$ тефо4 ю ргоуе а уапеу 
оЁ 1оБа[ сопуегоепсе гаёе 4ерепате оп Фе сопйригаНоп даа изе4 а{ еасВ Цегайоп $ер. 
Тре еззепсе оР1еасвтх соп5$155 ш сНоо$ше ай еасН Неганоп %ер засВ сопйгиганоп даа, 
УМсСЬ ргоу1ае а тахипит зрее4 оЁ сопуегоепсе. 


Статья поступила в редакцию 31.05.2012. 
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